ITEP-TH-45/12 



Bruhat Order in Full Symmetric Toda System 



12 



Yu.B. Chernyako 
G.I Sharygin 1,2 !^ 



, chernyakov@itep.ru 
, sharygin@itep.ru 



A.S. Sorin 2 , sorin@theor.jinr.ru 
Abstract 

In this paper we discuss some geometric and topologiacal properties 
of the full symmetric Toda system. We show by a direct inspection that 
the phase transition diagramm for the full symmetric Toda system in 
dimensions n = 3, 4 coincides with the Hasse diagramm of the Bruhat 
order of symmetric group 53 and 54. The method we use is based on 
the existence of a vast collection of invariant subvarieties of the Toda 
flow in orthogonal groups. We show how one can extend it to the case 
of general n. The resulting theorem identifies the set of singular points 
of dim = n Toda flow with the elements of the permutation group 
S n , so that points will be connected by a trajectory, if and only if the 
corresponding elements are Bruhat comparable. We also show that the 
dimension of the submanifolds, spanned by the trajectories, connecting 
two singular points, is equal to the length of the corresponding segment 
in the Hasse diagramm. This is equivalent to the fact, that the system 
at hand is in fact a Morse-Smale system. 
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1 Introduction 



Non-periodic Toda system (Toda chain) consists of n particles on the line 
with interactions between neighbours. The Hamiltonian of this system is 
given by 

n j n— 1 

H = 2^ + Yl exp( ^ ~ (!) 
i=i i=i 

where p, is the impulse of the i th particle and is its coordinate. The Poisson 
structure on the phase space (pi, qj) has the well-known form 

{p i ,qj} = 5 ij , {Pi,Pj} = } {qi,qj} = 0. (2) 

The evolution of the system is given by the standard Hamiltonian equations: 
p\ = {H, q[ = ~{H, Pl }. 

If we make the following ansatz ( j3]) 

bi=pi, ^ = exp-(ft - q i+ i), (3) 
the Hamiltonian will take the form 

n _. n— 1 

i=l i=l 

Observe, that there are only n — 1 variables aj. The Poisson structure (j2J) 
turns into 

aj_i} = -a»-i, aj} = a». (5) 

All the other brackets of coordinates and bj are equal to zero. In these 
coordinates it is easy to find the Lax representation of the system: one can 
show that the Hamilton equations are equivalent to the following matrix 
equation 

L=[L,A] (6) 
where L is called the Lax operator and is given by 
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and A is the operator 
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This system was first considered in [H [2], and in the work [3] there were 
found n functionally independent integrals of the motion. The involution of 
the integrals of motion was proved in the works [U EJ [7]. In particular, it is 
easy to see, that the trace of L is a Casimir function of the system, so that 
we can set it equal to 0, without any loss of generality. Thus we can assume, 
that Tr(L) = 0. 

There is a strightforward way to generalize these formulas: as one sees, 
the Lax matrix of this system is tri-diagonal traceless symmetric matrix, 
and A is the tri-diagonal anti-symmetric matrix, obtained from L by the 
natural procedure of deleting the principal diagonal and changing the signs 
of the elements above it. Now we can introduce a new system by taking L 
to be an arbitrary traceless symmetric matrix and A its anti-symmetrization 
and defining the dynamic by equation (Q. This equation is called the full 
symmetric Toda system. Below (see paragraph I2.1.ip we shall describe this 
system in more detail. 

As it follows from equation ($x§, eigenvalues Aj of L are invariants of the 
motion. Let us assume, that all these numbers are different. Moser in j6] 
showed that when t — > — oo the Lax operator of the usual tri-diagonal Toda 
lattice converges to the diagonal matrix with eigenvalues put in the increasing 
order, and when t — > +oo it converges to the diagonal matrix with decreasing 
eigenvalues. This property has been further studied in [12], where a general 
framework for calculating the eigenvalues of a symmetric matrix using the 
Toda equations was developed. The physical interpretation of this result is 
the full momentum transfer between n particles. The analogy of this property 
emerges in the system of (for example) two balls on a line. One of the balls 
move with the momentum p another ball has zero momentum. After elastic 
impact the balls change momentum. 

In |18| it was shown this property is true for the full symmetric Toda 
system and in [22] the transitions between the stable points of this equation 
in the case n = 3 for the full symmetric Toda system were described with 
the help of the explicit formulas for the solutions. 
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In our paper we study the behavior of the Lax operator of the full 
symmetric Toda system at t — » ±00. It turns out, that it always converges 
to a diagonal matrix with fixed eigenvalues, whose order at ±00 depends on 
the choice of the initial point (at t — 0) ot the trajectory. More accurately, 
we prove the following theorem: 

Theorem 13.11 The singular points of full symmetric Toda flow on SO n (M) 
can be identified with the elements of the symmetric group S n so that two 
points are connected by a trajectory if and only if they are comparable in 
Bruhat order on S n ( after suitable permutation, corresponding to the chosen 
order of eigenvalues). 

In this theorem we use the identification of the group S n with the 
group of permutations of the eigenvalues of L. Thus, we can describe all 
possible limits of the Lax operator: it turns out, that we cannot obtain 
any given permutation of the eigenvalues, when t — > ±00. In effect, the 
only permutations that are allowed are those, which go along the edges of 
the Hasse diagramm of the Bruhat order. Moreover, we can describe the 
dimension of the space of trajectories, that define such permutation: we prove 
the following 

Corollary 13.31 The dimension of the space, spanned by the trajectories, 
connecting two singular points is equal to the distance in the Hasse graph of 
the Bruhat order between these points. 

In fact, we show, that the system can be described as a Morse-Smale 
system on the flag manifold, so that its stable (resp. unstable) submanifolds 
are equal to the dual (respectively ordinary) Schubert cells. 

The rest of the paper is organized as follows: in section 2 we give synopsis 
of the facts from dynamical systems, Morse theory and the Schubert calculus, 
necessary for understanding the rest of the paper. In section 3, which consists 
of three subsections, we prove the main theorem. We begin with the simplest 
cases, n = 3 and n = 4 (paragraphs 3.1 and 3.2) and in 3.3 we prove the main 
results of this paper (theorem 13.11 and corollary I3.3p . Finally in appendix we 
describe the topology of the full flag space of R 3 . 

2 Definitions and conventions 

This section contains standard results and definitions from various 
mathematical theories, that we use in this paper. This section is also the 
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place, where we set notation and describe the conventions we use. 



2.1 Generalised Toda flows 

In these two subsections we shall recall the basics on the full symmetric Toda 
flow. 



2.1.1 General results 

As we have mentioned in the introduction, it is possible to generalize the Toda 
system (tri-diagonal Toda chain) and to consider the full symmetric Toda 
system. It turns out, that this system is in effect a Hamiltonian system: one 
can regard it as the dynamical system on the orbits of the coadjoint action 
of the Borel subgroup of SL n (M.) (equal to the group of upper triangular 
matrices with determinant equal to 1) see [El El EU fTTj . Namely, consider 
the following decomposition (we use the Killing form on sl n to obtain the 
identifications in the second and the third line): 

Sin = 50 n © b+, 

< = (KY © (son)* = Symm n © n+ (9) 
(b+)* = ( S o n ) x = Symm n , (so n )* - (b+)^ = n+ 

As one sees, we can identify the space of symmetric matrices with the dual 
space of Lie algebra of Borel subgroup: Symm n = (b+)*, and hence we can 
introduce a symplectic structure on Symm n , pulling it back from (b+)*. We 
define a Hamiltonian function on Symm n by the formula H(L) = Tr(L 2 ). 
Then the full symmetric Toda system is just the Hamiltonian system, 
associated with this data. 

For example, if n = 3 we shall have 

an 012 a i3 
L = | a 12 a 22 a 2 3 | then A 

y a 13 a 23 ^33/ 

where A = A(L) is the antisymmetric matrix, determined by L: A(L) = 
L + — L- (the upper half of L minus the lower half). 

It turns out, that the Hamiltonian system defined in this way is integrable. 
Some of the integrals of this system are easy to guess: take Hk(L) = 
Tr(L k ), k = l,...,n, then {Hi, HA = 0, which is proved by direct 
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computation, and H(L) = H 2 (L). However, this set of integrals is not 

fi(fi -)- 1) 

complete (in a generic point): dimension of our space is , and 

the number of integrals is only n; thus one should look for additional 
invariants. First they were found in [13] using the chopping procedure. 
Another (geometric) construction of the integrals was proposed by the first 
and the third authors of this paper, see [23]. This construction is based on 
the Plucker embedding of the flag space into products of projective spaces. 

In effect, similar constructions can be considered for arbitrary Cartan 
pairs (decompositions of a semisiple Lie algebra into a direct sum of its 
compact sublagebra and a subspace which verify certain properties). For any 
such decomposition one can associate an integrable system analogous to the 
Toda chain. In this paper we shall restrict our attention to the usual SL n 
case. A reader, interested in the general situation of arbitrary Lie group 
should refer, for example to the paper [17] . 

The evolution of L, determined by the Lax equation (see ([6])) can be 
described in the terms of orthogonal group: first of all, the equation suggests 
that the set of eigenvalues of L does not vary with time (although the order 
of these values need not be preserved) - this is equivalent to the invariance of 
traces L . Let A = diag(Xi, . . . , A„), ^2 Aj = be the matrix of eigenvalues 
of L. Then one can find an orthogonal matrix \1/ = ^(t), such that 

L{t) = tf^Atf-^t), tf(t) G SO(n,R). (10) 

The matrix \1/ evolves in time according to the law: 

(U) 

where now we put A = A(ty) is the composition of A(L) with the expression 

A(^f) = (\^Avp _1 ) + - (M/A^ _1 )_. (12) 

For instance, when n = 3, we have the following expressions for the matrix 
elements Oy of L and A in terms of the matrix elements ip^ of \I/ and the 
eigenvalues: 

n 

aij = Xk^Pik^Pjk- (13) 

k=l 
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It is also very convenient to write down the additional integral of the system 
in terms of ipij and Xf for n = 3 there is only one additional integral 

T . . ^13^33 . , , ^12^32 . . . V'llV'Sl 

J 1,1 — h A1A3 r A2A3 . 

013 013 a 13 

One can show, that this expression is in fact a Kasimir operator in this case. 

For any n the additional integrals can be chosen to be equal to the rational 
functions of appropriate minors of matrix Namely, they are equal 

l...k 

to the determinants of the submatrices of spanned by the intersections 
of first (or last) k rows and by arbitrary k columns of The equations of 
motion of these minors have the following form 

M' = f(\,^)M, (14) 

where / is some function of A and ip. 

So, the level set M = of M is an invariant subvariety of the Toda system. 
We shall sometimes call these surfaces the minor surfaces of the toda flow. 
For instance in the case n = 3 there are six such surfaces, which correspond 
to six 1 x 1 submatrices and when n = 4 there are eight distinct lxl minors 
and six 2 x 2 minors, which give invariant surfaces. Details can be found in 
|22] and |23] for the case of full symmetric matrix of the Lax operator, see 
also [13] where the Hessenberg matrices of the Lax operator are considered. 

2.1.2 Potential function 

An important property of the system discuss is, that on the level of the 
orthogonal matrices it is gradient, i.e. there is a function F(^>), such that 
its gradient, with respect to certain Riemannian structure on SO n (R), is 
equal to A§> . Observe also, that this function also depends on the matrix 
of eigenvalues A. In general case the corresponding functions and symmetric 
forms on so n are described, for instance, in [15|.|16] (for orthogonal group) 
and in [17] (for general case). It truns out that in the case we consider here 
the matrix A can be written in the form 

A=J([L,N}), 

where iV is a symmetric matrix, determined by the root system of SL n (M). It 
is not uniquely defined; for instance, one can take iV = diag(0, 1, . . . , n — 1), 
and J - an invertible symmetric linear operator on the space of antisymmetric 
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matrices so n , equal to the division by k on the fc-th upper- and lower- 
diagonals of a matrix. The operator J is later used in the definition of the 
Riemannian structure: one puts 

(A, B)j=(A, r\B)), 

where A, B G so n are arbitrary antisymmetric matrices and (A, B) is the 
Killing form on so n , i.e. in our case 

(A, B) J = -Tr(AJ- 1 (B)). 

Using this Riemannian structure we can represent A(^)^f as the gradient of 
the function 

F n (tf) = Tr(LN) = Tr^A^N). 

In fact, if we take ifr = (1 + 0)^ O) where G so n is an infinitesimal 
deformation of then we can linearize the function F n near \l/o : 

F n (*) = Tr((l + 0)^ o A^o '(1 " ®W 

= F(* ) + Tr(0^ o A^ o x iV) - Tr(^ A^o + °(@) 

= F(* ) + Tr(0[* o A*o \ ^]) + o(6) 

= F(* ) + Tr(ej-\ J([* A*o iV]))) + o(0) 

= F(* ) - (©, J([* A*o \ JV])> j + o(0), 

whence the result we need follows directly, since by a direct computation one 
sees that J[^ A^ \N} = A(*„), i-e. 0rad {)>J F n (tf o ) = -A(tf A*o ^o- 

It is instructive to write down explicit formulae for F n . For instance, if 
n = 3, we can take 

(-\ N 
N= 
\ 1. 

and 

F 3 (L) = Tr{LN) = a 33 - a n , 

(15) 

F 3 (A, ^) = Ax^ii - ^i) + A 2 (^ 2 2 - 1&) + A 3 (V4 - A 2 3 )- 
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2.2 Basics of the Morse theory 



In this section we recall the most basic definitions and results from the Morse 
theory. We warn the reader at once, that we shall rather need the geometric 
part of this theory, which deals with the local structure of the trajectories of 
a gradient vector field, and not the topological results, which constitute the 
major contribution of this theory to Mathematics. 

So let M be a compact smooth n-manifold without boundary. Recall, 
that a point x G M is called singular point of a function / : M — >■ R, if 
df(x) = (here df is the differential of /). Then / is called Morse function, 
if the following three conditions hold: 

1. / has only discrete set of singular points; 

2. the values of / in its singular points are different; 

3. all its singular points are non-degenerate. 

In effect, the first condition follows from the third, but we shall not need this. 
Recall, that a singular point x G M of / is called regular, or nondegenerate, 
if its Hesse matrix 

d 2 f d 2 f \ 

2 ' dxidx r 

d 2 .f 



( 



d 2 f(x ) = 



d 2 f 

dx'f 
d 2 f 



dx\dx2 
d 2 f 



dx2dx\ 



> d 2 f 

\dxndxi 



dx 2 , 



d 2 f 

dx n dx2 



dx2dx n 



dx 2 / 



\x=xo 



is nondegenerate in some coordinate system (and hence in any coordinate 
system) . 

Hesse matrix is a real symmetric matrix, and as one knows from linear 
algebra, its only invariant is the number of positive and negative squares in 
the canonic form of the corresponding quadratic function: 

n rj2 f n 



fe=i 



where 9k are new coordinates and Ek = ±1. It is this invariant, called the 
index of the point xo (which is equal to the number of —1 in the sequence 
(ei,...,£ n )) which plays crucial role in the whole Morse theory. 

An important stage in the development of the Morse theory is the 
following Morse lemma: 
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Figure 1: Types of singularities: a - source, b - saddle point, c - sink. 
a b c 

Lemma 2.1 Let f be a smooth function on a manifold M and Xq a 
nondegenerate singular point of f and k is its index. Suppose that there exists 
an open neighborhood U ofxo, which contains no other singular points. There 
exists a coordinate system (possibly on a smaller open neighborhood V oJxq), 
9\, . . . , 9 n , such that Xq has coordinates (0, . . . , 0) in this system and 



f\ v = f ( Xo )-el el + el +1 + ... + e 2 n . 



Let now g^ be a Riemannian structure on M. Then one can consider the 
gradient vector field of the function /: 

■ df 

dxj 

where g l i is the inverse matrix of g^. It is clear, that the singular points of 
this field coincide with those of the function (since the gradient is obtained 
by raising the indices of the differential). Moreover, the coordinate system 
8\, . . . , 8 n can be chosen so that the Riemanian structure gij is trivial in the 
point Xq (this result is sometimes referred to as the Morse-Palais lemma). 
Thus we can draw the local picture of the trajectories around a singular 
point (in fact the trajectories will correspond to the trajectories of the 
corresponding field on a Euclidean space up to an infinitesimal correction 
term). Thus in dimension 2 for instance we obtain one of the pictures from 
figure [TJ 

In general we see that the following statement holds 

Proposition 2.1 If the index (number of — l's in the canonic form) of a 
singular point is equal to k, then the trajectories that enter this point, span 
(locally) a k-dimensional submanifold. 
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An important particular case of the Morse system is the case of the 
so-called Morse-Smale systems. Namely, let us call the stable (respectivle, 
unstable) manifold of a gradient flow at a singular point p of this flow, the 
submanifold, spanned by all the trajectories, outgoing from (resp. incoming 
to) p. These submanifolds depend not only on the function, but also on 
the Riemannian structure and are usually denoted by W?(f,gij) (resp. 
Wp(f,gij)). Then we have the following definition 

Definition 2.1 One says that the gradient flow of a Morse function f 
determines a Morse-Smale system, if all the stable and unstable manifolds 
of this flow intersect transversely. 

In this case one can describe the dimensions of the intersections of 
the stable and unstable manifolds in the terms of the dimensions of these 
submanifolds and the dimension of the surrounding manifold M: 

dim(W;(f, gij ) f| W?(f, g l3 )) = dim W;(f, 9ij ) + dim W?{f, 9ij ) - dim M. 

This formula follows from the general results on transversality. This property 
makes the Morse-Smale systems important tools in algebraic topology. 

2.3 Bruhat order on permutation groups and Schubert 
cells 

In this section we shall recall the basics on the Bruhat order on the symmetric 
group and its relation to Schubert cells (for more details see, for example, 

m, to- 

Recall, that a partial order on a set X is a relation x < y on the elements 
x ytz y £ X, which verify the following condition: if x < y and y < z, then 
x < z. One of the convenient ways to visualize a partial order is by its Hasse 
diagramm: it is the oriented graph whose vertices are the elements of the set 
X, and the edges xy" correspond to the pairs x < y such that there is no z 
verifying the relation x < z < y. We shall say that the elements x and y for 
which this is true, are direct neighbours. Observe, that it is enough to know 
only direct neighbours to restore the whole partial order. The order, which 
appears in this way will be called the minimal order, generated by the given 
set of relations x < y. 

Let us now consider the symmetric group S n . It consists of all the 
permutations w : {1, 2, . . . , n} — > {1, 2, . . . , n}. We shall denote such 
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a permutation by the symbol w = (} ? ' ' ' U ) or simply by 

w = • • • 5 in), which means that w(k) = i^. One can abbreviate 

this notation to the so-called cyclic decomposition of w, i.e. one can 
represent an element w as a composition of non-intersecting cycles w = 
{ii,i 2 , ■ ■ ■,ik}{ji,32, where the symbol {n,i 2 , ■ ■ ■ ,ik} denotes the 

permutation 0, such, that <p{i\) = i 2 , 4>{i 2 ) — 13, ■ ■ ■ , (f>(ik-i) = ik, 4>{ik) = ii, 
and <j)(j) = j for all j 7^ i p , p — 1, . . . , k. In particular, (i,j) will denote the 
simplest 2-term cycle, which interchanges i and j. 

One defines the height (or length) of a permutation l(w) as the number 
of inversions in w: 

l{w) = #{(i,j) I 1 < i < j < n, w{i) > w(j)}. (16) 

One can now introduce a partial order on S n , the (weak) Bruhat order as the 
minimal partial order, generated by the following relations: 

x < y, if and only if y = (i,j)x, and l(y) = 1 + l(x). (17) 

Another approach to the definition of Bruhat order is based on the notion 
of Schubert cells or Bruhat cells and Schubert (Bruhat) varieties in the flag 
manifold. Recall, that a (full) flag E. in W 1 is a collection of hyperplanes 

{0}=£ C J E 1 CE 2 C---C E n _ x GE n = W n 

in R™, such that dimEj = i. The set of all flags constitutes a smooth compact 
manifold Fl n (R) of dimension (one can introduce coordinate systems 

on it from the evident identification Fl n (R) = SL n (R) /B+ = SO n (R) C\ B n)- 
For each permutation w, one defines the corresponding Schubert cell as 
the following subset X w G FL n (M): 

Definition 2.2 Schubert cell X w C Fl n (R) is the set of the flags verifying 
certain conditions: 

X w = {E. G Fl n I dim(E p r\F q ) = #{i | 1 < % < p, w(i) < q} V 1 < p, q < n}, 

(18) 

where flag E. = E l C E 2 C ... C E n = W l and dim(£ , i ) —i, F q — (ei, e g ) 
spanned by the first q elements of the basis ei, e n for M, n . We shall denote 
by r w (p, q) the number #{i \ 1 < i < p, w(i) < q}. 
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Schubert variety X w is denned as the closure of the corresponding Schubert 
cell, or which is equivalent by replacing all the equalities for dimensions 
by inequalities >. Then one can show that w < w' in Bruhat order, iff 
X w C X w >, see |19| (in the appendix, we describe a set of equations that 
determine the Schubert varieties in Fls(W), or rather the equations, that 
determine their preimages in (R)). 

In addition to the Schubert cells X w , one can define dual Schubert cells 
Q w , as the set of those flags E., for which for all p, q we have 

dim(E p p| F q ) = #{z < p I w(i) >n + l-q}, 

where F q is the subspace, spanned by the last q vectors of the basis. All 
the statements about the Schubert (or Bruhat) cells and varieties can be 
transferred to the dual cells and varieties with minimal modifications, see 
Fulton's book [19J, chapter 10 for details. Below we shall mention those 
modifications, that we shall need. 

There is a more general definition of Bruhat cells (and dual cells) and 
Bruhat order, which works for an arbitrary generalized flag space, i.e. the 
factor-space of a semi-simple group G by its maximal Borel subgroup B. 
Namely, one can choose an element w G G, representing some element of 
the Weyl group of G and consider the coset wB G G/B. Then Bruhat cell 
X w is equal to the orbit of wB under the left action of B in G/B. One can 
show, that this construction defines a cell decomposition of the flag space, in 
particular, one can show that these orbits are homeomorphic to open discs. 
The corresponding Bruhat order on Weyl group of G is then given by v < w 
iff X v G X w . 

In the case of the usual "geometric" flags, i.e. G = SL n (M.) the 
homeomorphism of X w and open disc of dimension l(w) (in fact, even IR^" 1 )) 
can be described in a rather straightforward way (see §10.2 of |19|): each 
flag E. has E p spanned by the first p vectors-rows of a unique "row echelon" 
matrix, the p th row of which has 1 in the w(p) th column with all 0's after this 
1, and the matrix has all 0's below these l's. For instance if w = (2143) such 
matrices have the form 

/* 1 0\ 
_ _ 10 
a ~ 0*1 ' 

\0 1 0/ 

where * denotes arbitrary real numbers. The corresponding Schubert variety, 
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however need not look like a submanifold, but one can describe this variety as 
the common zero set of a finite collection of polynomial equations on matrices 
in SL n (M) (or in SO n (M)), representing the flags. 

To this end consider a matrix A G SL n (M). As we mentioned before, this 
matrix determines a flag E. so that E p is spanned by the first p row-vectors 
of A (in some books column- vectors of A are considered). Let F. be the 
standard flag, i.e. flag, spanned by the rows of the unit matrix. In this case 
the dimension of the intersection E p f] F q is equal to the difference p — rk A Ptq , 
where A VA is the sub matrix spanned by the first p rows and the last n — q 
columns of A. On the other hand, p — r w (p,q) is equal to the rank of the 
similar submatrix of the permutation matrix W G O n (R), associated with w. 
Thus we obtain the following alternative definition of the Schubert cell: 

Definition 2.3 Schubert variety (closure of the Schubert cell) associated 
with a permutation w G S n is equal to the subset of flags in i<7+(R) = 
SL n (M)/Bn, spanned by all the matrices A G SL n (M.), such that all their 
north-eastern rectangular submatrices have ranks, less or equal to the rank 
of the corresponding submatrix of W G O n (M) where W is the orthogonal 
matrix associated with the permutation w. For dual Schuber cell one should 
compare the ranks of the north-western submatrices. 



Equivalence of these two definitions is proved in [19, 20J. 

(\ 2 3 4 

For instance, take permutation w = I ^ ^ ^ ^ 



then 



W 



/0 1 0\ 

10 

1 

\0 1 0/ 



and we can consider the matrix of ranks (i.e. the matrix, whose z, jth entry 
is equal to the rank of the rectangular submatrix Wij, spanned by the first 
i columns and j rows of W: 



rk(W) 



/0 
1 
1 

V 



1\ 

2 
3 



Now the condition that a matrix A defines a flag, that falls into the dual 
Schubert variety Q w is obtained by comparing the ranks of the submatrices 
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in A and the entries of this matrix, as we described above. Thus we should ask 
that the ranks of certain submatrices in A are less than what is prescribed by 
their sizes, which is equivalent to vanishing of all the minors, associated with 
these rectangular submatrices. For instance, in the case we consider here, we 
have the condition an = 0, det(a pg )i< Pi9 <3 = 0. 

In general, we obtain the polynomial conditions, saying, that all minor 
of prescribed sizes in a rectangular matrix vanish. Thus the Schubert variety 
is indeed a variety (on the level of matrices) and we have a complete set of 
equations, describing it (although it is, of course, excessive). 

We conclude this section with two general facts about the Schubert cells 
and dual Schubert cells 

1. Schubert cell X w and dual Schubert cell Q u intersect if and only if 
u < w in Bruhat order, in which case they intersect transver sally; 

2. Tangent space to X w at w is spanned by those positive roots of SL n (M), 
which are mapped into negative roots by w; tangent space to the dual 
Schubert cell Q w at w is spanned by those positive roots of SL n , which 
are mapped into positive roots by w (here we consider w as an element 
of the Weyl group of SL n (M) and identify the tangent space to the flag 
manifold SL n (M.)/B£ with the image of the tangent space to SL n (M) 
under the natural projection). 

Both this statements are well-known (see, for example, [19| . |21|). 

3 The phase portrait of the Toda flow 

In this section we shall give a description of the asymptotic behaviour of the 
phase trajectories of Toda system for low dimensions. The methods we use 
are rather simple-minded: knowing, that the system in gradient, in fact even 
Morse, we know, that every trajectory should flow from one singular point of 
the field to another (in particular, there are no periodic trajectories). Besides 
this, for each singular point of the vector field we can describe its local space 
of incoming and outgoing trajectories. Then we use a suitable set of invariant 
hypersurfaces, to describe the combinatorics of this transfer. 
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3.1 The case of S0 3 (R) 



We begin with the accurate description of the phase structure of the full 
symmetric Toda flow for n — 3. First of all we shall fix a matrix of eigenvalues 

Ai o o \ 

A = A 2 , Ai + A 2 + A 3 = 0. 
\o o X 3 J 

We shall assume, that Aj are all distinct, more precisely that 

A 2 > Ai > > A 3 . 

Choosing this order is equivalent to the choice of the positive Weil chamber 
As we know, the evolution of the system can be now described in the 
terms of the evolution of the corresponding system on SO^iM), see (TTTj) . It 
is easy to see, that the vector field A(ty) (see (IT2|) ) vanishes, if and only if 
the matrix \I/A\I/~ 1 is diagonal. Since the eigenvalues of A are all distinct 
and coincide with eigenvalues of ^A^f , and the matrix is orthogonal, 
we conclude that \I/A\I/ _1 can be diagonal, if and only if \I> is a matrix 
with positive determinant that permutes the vectors ±ei, ±e 2 , ±e3 where 
e ii e 2i 63 is the standard base of M. 3 . Such matrices form a finite subgroup 
5*3 of SOs(NL) (here we give a list of such matrices in S , 3 (R), symbols Si 
denote the cosets of the elements of S3 modulo the group N = §1 of diagonal 
matrices, see below): 



(\ o\ (-1 0\ 

si: 10, -10 

\o 1/ \ 01/ 

/ 1 o\ /o -1 0\ 

s 2 : -1 , 1 , 

\ 1/ \0 1/ 

/ 1\ /0 1\ 

s~ 3 : 10, 0-10, 

\-l 0/ \1 0/ 



-1 \ /l \ 

1 , 0-1 , 

0-1/ \0 -1/ 

1 \ / -1 \ 

1 0, -1 , 
00-1/ \0 0-1/ 

-1\ / -1\ 

1 , -1 , 

10 0/ \-l / 

(19) 
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-1 0\ I 0\ /l \ /-l 

•s, : I 1 I , [ 1 I , [ -1 I , [ 0-1 

0-10 

1 

.s, : I 1 I . | I | . | - 1 1 . | 0-1 

-10 

o o -r 

•s (j : I 1 I . | -1 | . | I | . |-10 

1 

here we denote by s, the classes of a finite subgroup S3 which correspond 
to the permutations to, of the eigenvalues Xj. There is a normal subgroup 
N = §1 C which consists of the diagonal orthogonal matrices with positive 
determinants, so that the factor-group S3/N is isomorphic to S3. As one sees, 
this subgroup is equal to the intersection of SO^M.) with the group of upper 
triangular matrices B£- As we have mentioned, this means, that the factor- 
space SOs(M.)/N is homeomorphic to the flag space Fl 3 (M). 

We are about to describe the asymptotic behavior of trajectories of the 
vector field A(^f) on SO^M.) (or on Fl 3 (R)), i.e. describe the trajectories, that 
go from one singular point of this field to another. To this end we shall use 
two important facts about this vector field. First, we recall that this field is 
the gradient of function Fsfif), see section l2.1.2l (we can regard this function 
as a function on Fls(M), since its value does not depend on the action of N). 
Thus, we can assume all the standard properties of the gradient vector fields, 
in particular, we at once see that there are no closed trajectories of this field 
and the value of the function F 3 always increases along a trajectory. Below 
we shall show, that all the singular points of the field are non-degenerate. 
Hence, we can use the results of Morse theory as well. 

Second, one of the important properties of the Toda flow is the fact, that 
there is a vast collection of invariant surfaces of this vector field, i.e. of the 
surfaces, that are tangent to the field A(^f). These surfaces are described in 
section 12.1.11 In our case there are 6 distinct surfaces, corresponding to 6 
entries of the first and the third row of which can vanish. We shall denote 
these surfaces just by ipij. Observe, that these equations are also invariant 
under the action of the group N. 

So, let us begin with the Morse theory. 

Proposition 3.1 All the singular points of the system are non-degenerate, 
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their index is equal to the height of the corresponding element in the Weil 
group (with respect to the chosen positive Weil chamber). 

Proof. This statement is proved in a general case in the paper [T7]. Here, 
the simplest way to show this is by a direct inspection of all 6 points 
§i, i = 1, . . . , 6, listed above (we can choose arbitrary representative in the 
corresponding coset for our purposes). 



Si : U * U , s 2 : * U U , s 3 




s 4 : u u * , s 5 : u u * , s 6 





where * takes the values ±1. 



It is convenient to choose in a neighbourhood of a point ^ G SOs(M) the 
coordinate system, induced from the Lie algebra SO3 in the neighbourhood 
of the unit and transferred to \1/ by right translation with Let 9%, 62, 63 
be the coordinates in SO3, i.e. let 





( 


01 


&3 




-61 









\-o 3 


-9 2 






be a generic element of SO3. Now in order to express the quadratic part 
of Fsfif) in the singular points Sj, i = 1,...,6 in coordinates, given by 
translation of (#1, 9 2 , #3), we replace the matrix f in a neighborhood of 
by 

# = ii + + ie 2 ^ + (e 2 ), (20) 

where by si we denote some matrix of the class si and o(0 2 ) stands for the 
sum of terms of degree 3 and higher in #j. Thus, we have 

L(«0 = ^a^ t = s iAsr 1 +[e,siA^-V©^r 1 e+^[e 2 ,s i As i - 1 ] + +o(e 2 ). 

(21) 



18 



Thus we can express by direct calculation the quadratic approximation of 
F 3 near Sj in coordinates 6 2 , 9 3 ). We have gathered the corresponding 
results in the following table: 



i 


F 3 (si) 


4^3 


Index 


Minors 


1 


A3 — Ai 


(Ai - X 2 )9\ + (A 2 - A 3 )^ + 2(Ai - A 3 )0§ 


1 


012,013, 031, ^32 


2 


A3 — A2 


(A 2 - \i)e{ + (Ai - \z)6l + 2(A 2 - A 3 )6»| 







3 


Ai — A3 


(A 3 - X 2 )9\ + (A 2 - X 1 )9 2 2 + 2(A 3 - Ai)0§ 


2 


■011, "012, V>32, ^33 


4 


A2 — Ai 


(A! - X 3 )9 2 + (A 3 - A 2 )0| + 2(Ax - A 2 )0§ 


2 


■012,-013,^31, ^33 


5 


Ai — A2 


(A 2 - X 3 )9\ + (A 3 - Ai)0! + 2(A 2 - Ax)^ 


1 


■011,-013,032,033 


6 


A2 — A3 


(A 3 - A 2 )0? + (Ai - A 2 )0| + 2(A 3 - A 2 )0| 


3 


011,012,031,033 



The second column of this table contains the values of F 3 at the given point, 
the third column consists of quadratic parts of the Taylor series at Sj. The 
nondegenracy of all singular points is now quite evident. Moreover, it follows 
from the choice of the order of Aj, that the values of F 3 in §i are all different 



F 3 (s~ 6 ) > F 3 (s 3 ) > F 3 (s 4 ) > F 3 (s~ 5 ) > F 3 (si) > F 3 (s~ 2 ). (22) 

So the function F 3 is Morse, the fourth column of the table contains the Morse 
index (the number of negative terms in the canonic form of the hessian) of 
the points. □ 
In the last column of this table we have given the list of all minor surfaces 
(see section [T2~j) . containing the corresponding singular point. As one knows, 
every trajectory of a gradient vector field connects (when the time tends 
to ±00) two different singular points. In addition, all the singular points 
being nondegenerate, we can use the Morse theory to describe (in terms 
of the local coordinates) the space, spanned by the incoming and outgoing 
trajectories (they correspond to the negative and positive eigenspaces of the 
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Hessian matrix). On the other hand, we know, that once in a minor surface, 
a trajectory cannot leave it. 

Let now T s ~ k ipij be the tangent space of the minor surface ipij at Sk 
and T^^ij its intersection with the positive /negative eigenspace of d^ k F 3 
respectively. We conclude, that the outgoing trajectories, falling into T^ ipij, 
should connect s k with another singular point in ip^, and the trajectories in 
Tf k ipij should come to sj. from a point in ipij. Thus, if we can cover the positive 
and negative eigenspaces inside the tangent space T s ~ k SOs(R) with T^ipij, 
we shall be able to show in what minor surface every incoming/outgoing 
trajectory of Sk lies. 

Besides this, we know, that every trajectory moves from lower values of F 3 
to higher values of this function. Now, suppose, that §i has only one outgoing 
trajectory inside ipij and the only singular point above it in ipij is Sk- Then 
we must conclude, that this only trajectory connects §i and Sk and we don't 
even need to know, if the surface ipij is singular or not. Combining these two 
simple observations, we can draw the full portrait of the system. 

We are going to complete this program. Let us begin with the following 
table, showing the incidence of singular points and minor surfaces. 



i 


ipu 


i>12 


1pl3 


1p31 


1p32 


^33 


S~6 


+ 


+ 




+ 




+ 


S3 


+ 


+ 






+ 


+ 


S~4 




+ 


+ 


+ 




+ 


S~5 


+ 




+ 




+ 


+ 


si 




+ 


+ 


+ 


+ 




S~2 


+ 




+ 


+ 


+ 





(23) 



The plus sign signifies, that the corresponding point belongs to the surface. 
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We can also describe the restrictions of the quadratic Hessian forms to the 
tangent space of the minor surfaces. To this end we should consider the 
linear relation on (0i, 62, #3) arising from the equation ipij = at a point Sj 
and restrict d 2 s .F^ to the corresponding subspace. The linear parts have the 
following forms: 



• si : 








S2 



• s~ 4 : 







01- 


* 










-01 • * 




• * 


-03 ■ * 







03 ■ * 


-9i 


• * 


02 ■ * 


-03 


• * 





0i- 


* 


03 ■ * 







02 ■ * 


-02 


• * 






03 ■ * 
02 ■ * 


1 ■ * 



-02 ■ * 



-0! • * 
-03 ■ * 



The equations of minor surfaces being simply ipij = 0, the corresponding 
linear relations take a particularly simple form 0j, = for certain k — 1, 2 or 
3, thus the restriction of d 2 g F 3 is given simply by removing this coordinate. So 
we sum up the results in the following tables (we replace the actual coefficients 
at the squares of $i by their signs for the sake of brevity) 



Si 


V'll 


4>12 


1pl3 


1p31 


4>32 


1p33 


S~6 




-ol el 




-el -el 




-61 61 


S3 


-el + el 


-el + el 






-el -el 


-61+61 


S~4 




+6* el 


02 Q2 

~ U 2 - U 3 


+6\ el 




+61-61 


S~5 


+ol el 




-61 + el 




+el el 


+61 + 61 


S~l 




+el + el 


-el + el 


-6\ + 6\ 


-el + el 




S~2 


+el + el 




+el + el 


+el + el 


+el + el 





(24) 
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This table allows not only find the indices of singular points on minor 
surfaces, but also describe the direction (in the terms of the local coordinates) 
in which the incoming/outgoing trajectories flow. In many cases it is enough 
to show, that two points are connected by a trajectory. For instance, we know, 
that there is only one outgoing trajectory from s~ 4 , directed along 6>i (see the 
first table). However, there are two potential endpoints of this trajectory, s~ 3 
or s 6 and we cannot choose one of them, if we don't use the minor surfaces. 
On the other hand, this trajectory goes along the surface ip^i, where there 
is no singular points between s 4 and s 6 , thus we conclude, that this unique 
outgoing trajectory of s 4 will tend to s~ 6 with t — > +00. 

However, this table is not enough to describe all the trajectories. In order 
to make the picture complete we should consider the intersections of minor 
surfaces. It is easy to describe them here: for instance, the intersection ip\2 H 
ipis (on the level of the group) consists of special orthogonal 3x3 matrices 
with zeros in the second and third places of the first row. Thus the second 
and the third columns of this matrix give a orthogonal basis of the subspace 
M 2 C M 3 , spanned by the second and the third coordinate vectors. The 
remaining (first) column should be equal to the unit vector, orthogonal to 
this subspace. Thus, we conclude, that the intersection on the level of the 
group consists of the matrices 



for some t G [0, 2ir]. Now passing to factor space modulo the action of 
diagonal matrices, we can fix the element in upper left corner equal to 
1 and factorize the remaining part by the group Z/2Z. This gives us a 
topological circle inside the factor space. Similarly, if we take the intersection 
like ipi2 H ipsi , we shall obtain a wedge of two circles (in particular, this means 
that the intersection is no more a submanifold) . Indeed consider the matrix 
which satisfies the conditions ip i2 = 0, ^31 = 0: 



This matrix defines an orthogonal basis in M 3 : j\ = {ipn, 1P21, 0}, f 2 = 
{0, ^22, ^32}, = {1P13, 1P23, ^33}- So, from the scalar product / 2 ) = 
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we see that there can be two distinct types of such matrices: 



^22 ^23 Or 1p 2 l lfj 23 
V>32 ^33 / \ ^32 ^33, 



And from /*) = 1 we get: 




\ /-l 
— sinti or I — costi sin^ 
costi / \ sintx costi 



in the first case and 




— sint 2 \ /— cost2 sini 2 X 
cos t 2 or I sint 2 cos £2 
-1 / \ 10 

in the second (here £1, £2 £ [0, 27r]). Each of these families determines a 
1-dimensional submanifold in SOs(NL) which is isomorphic to 1-dimensional 
group S02{M) = S 1 inside ^12 D "031- Thus, ipwOipsi is homeomorphic to 
the wedge of these two circles (on the level of Fl% (R) , where we can fix ±1 
to 1). 

It is easy to see that different critical points inside the double intersections 
might belong to different circles. For instance, in the example we considered 
above the first circle contains Si,at t\ = 0,7T and s 4 ,at t\ = ±7r/2, and the 
second circle contains §4, at t<i = 0, 27r and sq, at £2 = ±7r/2. 

In order to describe the phase transition from one point Sj to another point 
Sj inside the intersections of minor surfaces more visually, let us consider the 
following graph language. First of all recall that the infinitesimal behaviour 
of the gradient field of a Morse function at a singular point on a dim = 2 
surface (and hence on the minor surfaces we consider here) falls into one of 
the following three classes according to the index of the Hessian (see (pM)) ): 

• index is equal to (both coefficients at the squares of the Morse 
coordinates are +1), 

• index is equal to 1 (one coefficient is positive, another coefficient is 
negative), 



index is equal to 2 (both coefficients are negative). 
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In the case we consider here, Morse coordinates can be chosen proportional 
to {01,82,63) (up to infinitesimal terms), and hence all the the signs are 
determined by the signs of the coefficients at Of. We shall denote the 
corresponding vector fields symbolically by graphs with one vertex (singular 
point) and oriented edges, labeled with the coordinates 0j. The edge will be 
oriented towards the vertex if the corresponding coefficient is negative, and 
from the vertex in the other case. We shall also draw the lines, that leave the 
vertex above this vertex, and the incoming lines below the vertex. Thus we 
have the following three basic pictures 




Using these building blocks we obtain the following graphical analogue of 
the table §H§: 
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One can use this graphical language to describe the transfers in the 
intersections of the minor surfaces. Take for instance if) 12 PI ^13. From the 
pictures we see that this intersection contains only two points, §4 and s\. The 
next question is what trajectories are in the intersection of these two points. 
It follows from the graphical description of the surface ^13 = 0, that there is 
one trajectory in it entering s±, and that it goes along the coordinate 9\. Since 
the value of the Morse function in §1 is less than in §4, this trajectory does not 
belong to the intersection. Similarly, the trajectory, that leaves s\ along the 
coordinate #3 (see the picture for the surface ipi2 = 0) does not belong to the 
intersection, since it does not belong to the surface from t/> 13 = 0. Analogously 
for §4 we must delete the incoming trajectory along the coordinate #3 and 
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the outcoming trajectory along the coordinate 9\. 

The following table sums up these considerations, it describes all possible 
transfers from one point Si to another point Sj which appear in the 
intersections of various minor surfaces ipij fl ipkf- 



3 


V>n 


0ia 


■013 


031 


V"32 


■033 


011 


X 


s 3 -> s 6 


S 2 -4 S 5 


S 2 5 6 


S 2 — >■ S5 , S5 — > S3 


S5 — > S3, S3 — > Sq 


012 


X 


X 


Si -► s 4 


Si — >■ S4, S4 — >■ Sg 


5i -> S 3 


S4 — > Sg, S3 — > Sq 


013 


X 


X 


X 


S 2 — y Si j Si — ► S4 


S2 — *■ Si, S2 — > S5 


S5 — > S4 


031 


X 


X 


X 


X 


S 2 Si 


S4 — *■ Sq 


032 


X 


X 


X 


X 


X 


S 5 -¥ S3 


"033 


X 


X 


X 


X 


X 


X 
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Finally observe, that however complicated the geometry of minor surfaces 
be, they are 2-dimensional surfaces in generic points. Thus when we restrict 
our attention to a surface of this sort, we see that a generic trajectory on it 
should flow from the lowest to the highest singular point lying in this surface, 
there being only a one-dimensional subspace, spanned by other trajectories. 
Similarly, a generic trajectory on SOs(M) / N should flow from the minimum 
to the maximum value of F3. Summing up, we conclude that the trajectories 
connect the singular points as described in figure [2j Here on the following 
figure thin black arrows denote the "most singular" trajectories, "fat" blue 
arrows correspond to two-dimensional flows, "the fattest" red arrow is the 
generic trajectory on SOs(B.), i.e. such trajectories span a three-dimensional 
subspace. As on can readily see, the structure of the arrows reproduces the 
Hesse diagramm of the group S3, the Weyl group of SL 3 (R). Moreover, the 
black arrows correspond to the minimal edges of the Hasse diagram, the 
blue arrows connect the points separated by one intermediate point in this 
diagram, and the red arrow connects the points, between which there are two 
other points. 
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Figure 2: The trajectories in SO^iM) 



3.2 Trajectories on SO^{ 



We now describe the trajectories of Toda system on SO^R) (or rather on 
S0 4 (R)/S0 4 (R)f)B+ ^ Fl A (R)). The ideas we use here are the same as 
before, but the quantity of calculations is much bigger (in particular, there 
are 24 singular points, which can be identified with the elements of 5*4). So, 
we begin with fixing the order of eigenvalues: 



A 2 > Ai > A 3 > > A 



The Morse function F 4 is equal to F±(a) = Tr(LN), where 



N 



/o 0\ 

10 

2 

\0 3/ 



Once again one can show that the singular points of this function on SO^{ 
coincides with the group of 24 • 8 = 192 matrices with det = 1, representing 
permutations of the vectors ±ei, ±e2, ±e3, ±&4. This group contains normal 
subgroup iV = A/4 of diagonal matrices (maximal torus), so that the phase 
space of our system can be identified with SO^(M.)/N, and the singular 



points of the induced system consists of 24 points 



where i^isH 



is a permutation of four elements, listed below (in orthogonal group * should 
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replaced with ±1, so that the determinant be equal to 1): 
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F, Morse 

Figure 3: Singular points and their indices in SO^( 
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As before, we introduce the local coordinates at a point g G SO^yK), pulling 
them from the Lie algebra so 4, whose generic element 6 is equal to 







/ e 1 

-9 l 



V 



-<7 3 
-04 



-<7 2 

-e 5 



o 3 e 4 \ 

9 % 

-e 6 ) 
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We can now express the second order approximation of F 4 at a point 
in the terms of coordinates (#1, . . . , 9 e ), for instance: 

^ 4 (5i234) = e\ (Ai - A 2 ) + 29\ (Ai - A 3 ) + 9\ 

+ (A 2 - A 3 ) + 39j (Ai - A 4 ) + 29 2 h (A 2 - A 4 ) + 9\ (A 3 - A 4 ) . 
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1234 
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Figure 4: Examples of £03 (R) clusters inside £04 (R) 



To get the corresponding approximations for the rest of s^i^i^, one should 
permute the eigenvalues A, by Thus we get the following picture of 

the singular values and their indices distribution (see figure[3]). At this picture 
we have shown all 24 singular points, vertical axis shows the indices of these 
points and the horizontal axis - the corresponding values of the function F4. 

As one can see, there are four clusters of six elements each, the points 
at every cluster resembling the points in the phase portrait of £03 (R). 
These clusters do, actually correspond to the subspaces, isomorphic to 
£03 (R) within £04 (R): the lower left cluster here is a genuine subgroup 
inside £04 (R), isomorphic to £03 (R), the subgroup, consisting of orthogonal 
transforms which leave invariant the fourth basis vector. There are four such 
subgroups (one for every vector e,, % = 1, . . . , 4). We shall denote these groups 
by £Og(R), % — 1, . . . , 4. It is straightforward to see, that these subgroups 
are invariant with respect to the dynamic system under consideration. For 
instance, one can prove this by the direct calculation with formula (ITOl) . 
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3412 



4312 

1(3)-3(2) > 
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4132 
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3142 1432 3421 4231 4123 



3124 3241 1342 2431 



1423 



4213 
" 4 



+2 



3214* 1324 



1243 



2341 



2413 



+ 1 



2314 



1234 



2143 



2134 



Figure 5: An exotic SC^R) cluster inside 5O4 



Moreover, this calculation shows, that the dynamic system, induced on these 
submanifolds from SO^R) coincides with the full symmetric Toda chain on 
SOa(fil), if we use the natural identifications (just observe, that the condition 
^2 Xi = is not necessary for the definition of the field A, see equation (flO]) ). 
Thus, the restriction of the flow diagram to these subgroups coincides with 
the diagrams of Bruhat order for SC^R). 

These observation can be further generalized to get the full description 
of the phase flow on SO^iM.). Let s be one of the matrices, listed above (see 
p. 27). Suppose that the permutation, corresponding to it, sends 1 to i\. We 
can assume, that the lii-th matrix element of s is equal to 1, and not to — 1. 
Then consider the coset space 

SO 1 / 1 =SOl{R)-s. 

This set consists of all the orthogonal linear operators,which send e\ to e^. 
It is evidently diffeomorphic to the group SOs(R), although it does not have 
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a natural group structure. 

Let now \1/ E SOl n be an arbitrary element. As it follows from the 
equation ffTUj) . the vector field A = A(^,A) at this point is equal to the 
right translation by s^j^^ of the field A(^s~ l , sAs" 1 ) on SO$(M), which 
can be identified with the vector field A on 5*03 (R), provided, we take the 
eigenvalues equal to a permutation of A2, A3, A4. Thus, the phase structure of 
the restriction of the Toda system to this subspace is equal to the translation 
of the corresponding phase portrait from SO3 (R) via right multiplication by 

The same construction is applicable if we replace 1 with 2, 3 or 4. Thus 
we get a collection of 16 intersecting clusters inside the set of all critical 
points of the Toda system for S , 4 (R), such that the restrictions of phase 
portrait to any of these clusters coincide with the portrait of 5*03 (R) (after 
suitable identifications). In table (|27|) below we give a list of these clusters. 
And picture H] shows the corresponding phase pictures: one, corresponding 
to S0\ A and the other given by SOf = SOl W- Observe, however, that the 
indices of singular points inside a subspace SO3 and inside SO^R) need 
not coinside. For example, at picture we have drawn the scheme of the 
cluster SOf-. As one can see, three of the arrows, which correspond to the 
singular points whose Morse index on the level of .SC^R) differ by 1, at this 
picture connect points, whose Morse indices in .SO^R) differ by 2. Moreover, 
the same pair of singular points can represent neighbors in one cluster and 
non-neighbors in another. 

Summing this up, we get the following table showing all the .SC^R) 
subspaces inside SO^R): here we use the notation, introduced above to 
describe the subspace we deal with; further, in the second column we give the 
list of singular points, that fall in this subspace; finally, in the last column 
we list several minors conditions, that single out the subspace. 

As one can show by direct inspection, all the dimension lesser than 4 
eigenspaces of Hessians of all the singular points in .SO^R) are covered by 
one or another of S0% subspaces listed in table (127|) . Thus we only need 
to check the configurations of minors surfaces with co-dimensions 2 or less. 
This is done by a direct checking. Thus, we obtain the following picture of the 
flows inside SO^M.) (or, rather, inside the corresponding flag space Fl^M.)), 
see figure |6j Here we show only the one-dimensional flows on 5*04 (R), that 
is the flows, corresponding to 1-dimensional eigen-subspaces of Hessians at 
singular points (for references we have added the cyclic decompositions of 
the corresponding permutations). 
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Cluster 


Points 


Minors 


sof 


4312, 4132, 4321, 4123, 4231, 4213 


1pll,1pl2,1pl3,1pM 


sof 


3412, 3421, 1432, 1423, 2431, 2413 


014, 044, Ml2 . Ml2 , Ml2 
12 13 23 


sof 


3142, 1342, 3241, 2341, 1243, 2143 


014, 044) Ml2 , Ml2 , Ml2 
14 24 34 


sof 


3124, 1324, 3214, 1234, 2314, 2134 


014,041,042,043 


sof 


2431, 2413, 2341, 2143, 2314, 2134 


011, 013, 0^14,042 


sof 


4231, 4213, 3241, 1243, 3214, 1234 


012, 042, Ml2 , Ml2 , Ml2 
13 14 34 


sof 


4321, 4123, 3421, 1423, 3124, 1324 


012, 042, Ml2 , Ml2 , Ml2 
12 23 24 


sof 


4312, 4132, 3412, 1432, 3142, 1342 


012,041,^43,044 


sof 


3412, 3421, 3142, 3241, 3124, 3214 


011,012,^14,043 


sof 


4312, 4321, 1342, 2341, 1324, 2314 


013, 043, M12 , Ml2 , Ml2 

12 14 24 


sof 


4132, 4231, 1432, 2431, 1234, 2134 


01 3 ,043,Ml2,Ml2,Ml2 
13 23 34 


sof 


4123, 4213, 1423, 2413, 1243, 2143 


013, 041, ^42, 044 


sof 


1432, 1423, 1342, 1243, 1324, 1234 


012,013,^14,041 


sof 


4132, 4123, 3142, 2143, 3124, 2134 


0H,041,Ml2,Ml2, Ml2 
23 24 34 


sof 


4312, 4213, 3412, 2413, 3214, 2314 


01 1 , 04 1 , M 12 , M 12 , M 12 


sof 


4321, 4231, 3421, 2431, 3241, 2341 


011,042,^43,044 
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Figure 6: The 1-dimensional flows at S'O^R) 



As one can see, it is just the Hasse diagramm of the Bruhat order on the 
symmetric group 5*4, or rather its image under the action of the permutation 
(\ 2 3 4\ 

1 2 1 3 4 j ' ^ S * S ^ 6 P ermu t a ti° n ' which reorders the eigenvalues Aj so 
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that Ai < A2 < A3 < A4. Thus we get for SO^iM), just as we did for .SC^IR) 
before, again the following result: 

Proposition 3.2 Two singular points x and y of the full symmetric Toda 
flow on SO±(Jxl) are connected by a trajectory if and only if x < y in Bruhat 
order of S4 (after suitable permutation, corresponding to the chosen order of 
eigenvalues). Moreover, for any n = 2,3,4,5,6 the n- dimensional flows of 
the Toda system correspond to the length n segments of Bruhat order (recall, 
that a segment in partial ordered set is a sub poset [x;y] = {z G P\x < 
z < y}, it is said to have length n, if for any chain of the maximal length 
x = Xq < X\ < ■ ■ ■ < Xk = y inside it, k = n). 

3.3 The general case 

Let us now prove that in the case of arbitrary n the similar conclusions should 
hold, namely: 

Theorem 3.1 The singular points of full symmetric Toda flow on SO n (M) 
can be identified with the elements of the symmetric group S n so that two 
points are connected by a trajectory if and only if they are comparable in 
Bruhat order on S n ( after suitable permutation, corresponding to the chosen 
order of eigenvalues). 

Moreover, we shall also prove that there is a fc-parametric family of such 
trajectories, if one can insert k points in Bruhat order between the given two 
(see corollary 13. 3p . 

The remark, concerning the order of the eigenvalues means that if we 
identify S n with the group of permutations of the eigenvalues of L so that 
the standard (default) order of the eigenvalues is Ai < A2 < • ■ • < A n , then 
the corresponding phase portrait can be identified with the Hasse diagramm 
of the Bruhat order. Otherwise, one should multiply the elements of S n in 
Hasse diagramm by the corresponding permutation. In previous sections we 
have described the low- dimensional case (n = 3, 4), where this statement can 
be proved by direct inspection. It is our purpose now to prove this statement 
for arbitrary n. 

3.3.1 Supporting remarks 

Before we give the full proof of the theorem 13. 1[ let us give few observations, 
supporting it. We begin with the following fact: if an element a G S n , n > 3 



35 



precedes (5 with respect to Bruhat order, so that there is no other elements 
of S n between them, then there is an element i G {l,...,n} such, that 
a(i) = /3(i). This is an almost direct consequence of the definition of Bruhat 
order (see section [2T3|) . since such neighbors always differ by a single swap of 
two elements. 

Now we reason by induction on n: suppose that for all k < n + 1 we know 
that any two elements in that are direct neighbors in Bruhat order (i.e. 
there are now other elements between them) are connected by a trajectory 
of Toda flow inside C SO k (R)/SO k (R) f| B+ = Fl k (R). 

Consider the singular points §j G S n+ i of the system in dimension n + 1. 
We shall always regards the elements of S n+ i C FL n (M) as the classes of 
equivalences of matrices inside SO n (R); similarly, we shall speak rather about 
the submanifolds in SO n+ i(R), and not about their images in Fl n+ i(M). 
Thus, one can consider (n+ l) 2 submanifolds SO% C SO n +i(M), isomorphic 
to the group SO n (M): they are defined in the same way we have previously 
defined the submanifolds SO$ C SO^(M.)/ SO^(R) f]B^. In particular, there 
are n + 1 subgroups SO % n {R) C SO n +i(R), isomorphic to SO n (R). As before, 
it is evident, that the restriction of the Toda system to these subgroups can 
be identified with the Toda system on SO n (M), and hence, the trajectories 
of the restriction of Toda system to any of SO l n (R) connect two points inside 
the corresponding part of S n if and only if they are comparable with respect 
to the Bruhat order in S n . The same is true for the restrictions of this system 
to all the other submanifolds SO?* C SO n+ x(R). 

On the other hand, by definition, we have the following characterization: 
any two elements a, (3 G S n +i belong to one submanifold SO l J for some 
i,j = 1, . . . , n + 1 if and only if there exist an element 1 < i < n + 1, such 
that a(i) = j = f3(i). Hence, every two points a, (3, neighboring with respect 
to the Bruhat order in S n+ i, will fall to the same SO% for certain Since 
Bruhat order is hereditary in the sense that the restriction of Bruhat order 
on S n C S n +i (which is embedded as the set of all permutations, preserving 
certain element) is induced from this order in S n+ i, we conclude, from the 
observation above, that they shall be connected by a trajectory of Toda flow 
onSO n+1 (R). 

Thus we have proved, that any two neighbors in Bruhat order will be 
always connected by a Toda trajectory. However, we still do not know, what 
happens, if these two points are not neighbors in this order. And we cannot 
tell anything about the dimension of the space of all trajectories, connecting 
any two points. 
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One can try to use similar ideas to prove the reverse statement (i.e. to 
show that the points, that are not comparable with respect to the Bruhat 
order, cannot be connected by a trajectory). Instead of this, let us proceed 
by proving the general statement. 

3.3.2 Proof of the theorem I37T1 

Proof. As we have mentioned, Toda system preserves varieties, given by the 
equations M i...k = 0, where M i...k is the determinant of the submatrix of 



spanned by the intersections of the first k rows and the columns ii, . . . , i^. 
Hence, if equations, determining the Bruhat cell are of this form, then this 
Bruhat cell is preserved by the Toda flow, and any point in such cell B a can 
flow only to a point, that lies in the closure of this cell, which consists of the 
cells, corresponding to those elements of S n +i comparable with a in Bruhat 
order. 

However not all cells are spanned by equations of the necessary form, since 
there can be rectangular subalgebras of type k x I with k < I (see section 
12. 3p . Thus we shall need the following proposition to be able to handle all 
Bruhat cells: 

Lemma 3.1 If the (upper-left) submatrix ty l,J (0) o/\l/(0) has rank less, than 
min(i, j), then the same is true for all t: rk \E rij (t) < min(i, j). 

Proof. The proof is based on the following observations: first, to say, that 
rk^'-^t) < min(z, j) is equivalent to impose the condition that all the minors 
of the min(i, j) x min(z, j) submatrices in ^(t) are equal to 0. If i = j or % > j 
(i.e. the matrix is square or "horizontal"), then the conclusion of the theorem 
follows at once from the fact, that equations M i...k (t) = hold identically, 



if they are verified for t = 0. If i < j, we need to work a little more. 

Namely consider the minors of which are equal to the intersection of 
the first k < n columns of \1/ with an arbitrary set of k rows. Then a direct 
computation shows that the following equations hold 



k k 



1,2,...,? 



k 



1,2 k 



1=1 j = l 



(28) 



k i<i s 



2 EE°^ (-ir (p) ^i 



1,2,...,*! 



s=l i=l 
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where ii...i...i s ...jk is a permutation p of the set ii, i, ...i s , jk, &{p) is 
the parity of p, i.e. the number of inversions in p. The indices in the sets 
ix, 12, ik and ii, i, ...i s , are increasing (recall, that the numerator a 
in the index | corresponds to the rows and b to the columns of and aj r i s 
is an entry of the Lax matrix L). To prove this formula consider the minor 
M i 1 ,i 2 ,...,i k . It is equal to the following polynomial 

1,2,...,*; 

M H , 2 ,..., k = E^-l^tyi^Vw-Vw (29) 

1,2,...,* y y ' 

where Jij2---Jfe is a permutation p of the set 1, 2, and a(p) is the number 
of inversions in this permutation. Differentiate this expression with respect 
to t using the formula 

k=i 

ip'ij = (-an + - 2 ^ a-k-u^k-ij- (30) 

k=l 

Change the subscripts so that the formula for dynamics ip takes the following 
form: 

i<i a 

^isj s = (~ a isis + Aj-JVi.i. - 2 (si) 

1=1 

Finally take the formula 

M ' n ,, 2 H , = E(- 1 )* (P) «ii^-^ 



i-kjk 

1,2,. ..,fc 



(32) 



+ VWiVw--^* + ••• + ^hh^h-AuJ' 

substitute (}3~Tj) in (132]) and after regrouping we get 

k k 

M'i^,...^ = ^(-l)' (p) ((~ J2a kil +J2x j )^ hjl ij; i2h ...ip ikjk 
i,2,...,* — 7-7 r - : 

P 2=1 J=l 

«<«i i<«2 



+ (-2 y ^2a iil il) ih )il) i2j2 ..4i kjk + i) iljl (-2^2a ii2 il) i j 2 )il) i3h ..4, 

i=l i=l 
+ ••■ + V ) njiV ) i2j2-"'0ifc-lifc-i(~ 2 ^ ^ a Hk' l l J ijk))' 



i=l 

(33) 
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which is equivalent to (j28p . 

Now formula ( 12 8 p shows that the derivative of M i 1 ,x 2 ,..., lk can be expressed 

1,2,. ..,k 

in the terms of the minors Mj j ^ , where at least one of the entries i', is 

1,2 fc 

less then the corresponding i\. Hence we can use induction to show that the 
equations M(t) i 1} i 2 ,...,i k = hold identically for all t. □ 

1,2 k 

The following statement is a direct corollary of this lemma: 

Corollary 3.1 Toda flow preserves Schubert cells (in other words, the 
corresponding vector field is tangent to Schubert cells). 

In fact, the lemma shows, that the system of equations, which define the 
Schubert cell is preserved by the Toda flow: tbhis is evident for the case, 
when the rank of the "vertical" submatrix differs from the maximal by 1. In 
a more general case, we should apply this lemma to the matrix, equal to \l/ 
with few omitted columns. 

If we define the dual Schubert cells as the sets of matrices, for which 
the conditions, similar to that of the definition 12.21 holds for the lower left 
submatrices, (see section I2.3P then we obtain a similar statement 

Corollary 3.2 Toda flow preserves dual Schubert cells. 

Now we can use the general facts about the Schubert cells and dual Schubert 
cells, listed in section 12.31 

We shall also use the following general proposition: 

Proposition 3.3 Let (M, g) be a riemannian manifold, N be a manifold, 
immerged in M via a smooth full rank map F. Let f be a function on M. 
Let g' and f be the restrictions of g and f to N . Then 

D F (grad g ,f') = proj^(N)(grad g f). 

Proof. The statement of this proposition is local, it is enough to prove it for 
the tangent space at a point, so we can assume, that M = M. a , N = M. b and F 
is a linear embedding. Choose orthonormal base for g' in R 6 and complete it 
to an orthonormal basis in R a (with respect to g). Then in the corresponding 
coordinates grad g f is just the vector of partial derivatives (J^-, • • • , ^r) and 
Dp(grad g i f) = (J^r, • • • , J^, . . . , 0), so that the statement is evident. □ 
We shall apply this statement to Schubert cells. One more statement, 
that we shall need is the following: 



39 



Proposition 3.4 If the order of the eigenvalues in A is given by X\ < A2 < 
• • • < A n; then the negative (resp. positive) eigenspace of the Hessian of the 
Morse function F n of at w is spanned by those positive roots of SL n , which 
are mapped into negative (resp. positive) roots by w. 

This proposition (in a general case of arbitrary Cartan pair) was proven in 
the paper [T7]. For the sake of completeness we give here our prove of this 
fact in the case of the pair Sl n (M) D 5 , O n (M). 

Proof. First of all let us fix some matrix s k G S n (the definition of S n is 
similar to the definitions of S3 and S4) 



( * 



* 

\ 





* 



\ 




/ 



(34) 



The non-zero elements * of Sk are equal to ±1. Conjugation with Sk defines 
a permutation w on the eigenvalues in the diagonal matrix A. 

The matrix \l/ in a neighborhood of s& has the following decomposition 



s~k + Qs k + -Q 2 s~ k + o(6 2 ) 



(35) 



Here o(6 2 ) denotes the sum of terms of degree 3 and higher in 9{j. Thus, we 
have 



L{9) = *A^ 3 



1 



s~ fe As fe s~ k As~ k ~ 1 }-Qs k As~ k - 1 Q+- [0 2 , s k As~ 



V 1 ] 



where the matrix O G so(n) has the following form 



(36) 



/ 

—#12 



e 



V 



"01n-l 



6\2 


—®2n 





~@n-ln 



@2n 




(37) 



We shall be interested in 



-&s k As~ k 1 e + -[e 2 ,sfcAs- fc x ] 
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to calculate the canonic form of the Hessian d 2 ~ k F n . Consider the diagonal 
matrix A w 

SkAsT 1 = K (38) 

then its matrix elements are Xa = A TO (j). We need only diagonal parts of 
the terms — QskAs^ 1 ® and + |[0 2 , SfcAifc ]+ because F n = Tr(LN), where 
diagonal matrix N has the matrix elements Nu = i — 1. Put 



(p x = - Qs k As k 1 6, (fl = 



3+i 

for j = 1, . . . , n. Here we replace 0y by 9ji if % > j. We compute 
dlF n = Tr((^ + cp 2 )N) 



i<j 



+ E©^)-^))(i-i), (39) 

i<j 



i<j 



Thus we have expressed by direct calculation the quadratic approximation of 
F n near s k in coordinates 6^ . The statement of the lemma now follows from 
a direct inspection of this formula. □ 

Finally, let us put together all the statements: it follows from the 
propositions and 13.41 and the second property of the Schubert cells, that 
Schubert cell X u (resp. the dual Schubert cell Q u ) coincide with the unstable 
(resp. stable) submanifolds of the Morse function F n at w. Indeed, from 
pr op osit ion 13 .3 .21 it follows that the Morse flow of the restriction F u \b u (resp. 
F n \D u ) coincides with the Toda flow. And from proposition 13. 41 it follows, that 
w is local maximum (resp. minimum) on B w (resp. D w ). Now the proposition 
follows from the first property of the Schubert cells and their duals. □ 

Finally observe, that we have actually showed, that Toda system in this 
case is a Morse-Smale system (see the definition 12. ip . Thus we obtain one 
more statement: 
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Corollary 3.3 The dimension of the space, spanned by the trajectories, 
connecting two singular points is equal to the distance in the Hasse graph 
of the Bruhat order between these points. 

Proof. The statement follows from the transversality of the Schubert and 
dual Schubert cells and the fact, that the dimension of the Schubert (resp. 
dual Schubert) cell of w is equal to the number of the elements, connecting 
it to the bottom (resp. top) element in the Bruhat diagramm. □ 

Remark 3.1 Another approach to the question of describing the phase 
transfer of Toda flow was given in the paper [22], in which the authors use 
explicit formulas for the solutions to describe their asymptotic behaviour. 
One of the outcomes of their approach is a construction of a polytope, whose 
vertices (after suitable identifications) correspond to the elements of the Weyl 
group. One may ask, if the edges of this polytope reproduce the structure of 
the minimal edges of Bruhat order (one can see, that this is so for SO3 (M) 
case, which is proved in the cited paper). 

Appendix 

A Geometry of the Toda flow on SO (3) 

In this additional section we try to visualize the results, concerning the lowest 
dimensional case we considered here, that of 50(3). 

This group contains a normal subgroup N of diagonal matrices (the first 
row of this list), which commutes with the matrix A. It follows, that the 
vector field M( 1 i r ) does not depend on the choice of element in the right N- 
coset of SOs(M). So we can consider the system as an evolution of points on 
factor-space SO^M.) / N . The singular points of this system in this setting will 
consist of 6 cosets (s i( i = 1, . . . ,6), which can be regarded as the elements 
of the permutation group S3 = W(sOs). one can visualize the structure of 
this space and positions of the singular points on it, using the quaternionic 
presentation of SO^Wi): 

S0 3 (R) = Hi/(Z/2Z) 

Here Hi denotes the unit quaternionic sphere and Z /2Z acts on it by ±1. 
Then the factor space SO^(M.)/N is equal to the factorization of S 3 by the 
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Figure 7: The Voronoy set of N 



preimage N of N in HI. One can find the explicit quaternions, representing 
the elements of N: 

±1, ±i, ±j, ±k, 

then the factor-space Mi/N can be identified with the factorisation of the 
spherical Voronoy set of 1 G HI (with respect to the other elements of N) 
along its boundary, see figured As one sees, this domain is (combinatorially) 
a cube. The action of N identifies the opposite faces of this cube after having 
turned them 90° around their centers. It is also easy to describe the positions 
of the singular points of M at this set, see figure [HJ at this figure there are 
(after identifications) shown all 6 singular points: 1 point is in the center of 
the cube, 3 points are the centers of the three pairs of opposite faces of the 
cube, and 2 more points are given by the orbits of the cube's vertices under 
the action, which identifies the opposite faces, turned by 90°. 

There are six Schubert cells X w in this case and we can write down the 
systems of equations, that determine each of them. In fact, we have a cellular 
decomposition 

FL 3 (R) = |J X w (40) 

of the flag manifold. As we have explained in section 12. 3[ we can describe the 
preimages of these cells in SL 3 (R) with the help of systems of equations on 
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Figure 8: The fundamental domain of N, positions of the singular points 



matrix entries (these preimages are often called the matrix Scubert varieties). 
Here we describe the intersections of these preimages with SO^M.) C SX 3 (1R). 
For each permutation io, we have the orthogonal matrix sf 



si : U * U , s 2 : * U U , s 3 






s 4 : u u * , s 5 : u u * , s 6 



where * takes the values ±1. The corresponding rank matrices and 
polynomials cutting out the matrix Schubert varieties inside £703 (R) have 
the following forms: 



rk(sx) = I 1 2 2 , no polynomials, 

\l 2 3/ 

/O 1 1\ 

rA;(s 2 ) =1 2 2 , = 0, 
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rk(s 5 ) 



rk(s 3 ) 



rk(s e ) 




l/j U = 0, ^12 



Ipll = 0, 1p21 



i>w = o, V12 



M12 = 0, 



12 



0, Mia 



0, V21 = 0, M12 = 



0, M12 = 



12 



12 



0. 



As one can see, in many cases the recipe, we use, produces more equations, 
than one really needs to describe the cell. Thus the dimensions of the the 
cells should be calculated with care: this table does not give them "on the 
nose". 
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